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ABSTRACT 

Quantized  Frequency  Modulation  is  one  of  the  various 
methods  used  to  reduce  the  effects  of  multipath  propaga- 
tion of  high  frequency  radio  waves .   The  modulation 
technique  is  described, and  a  signal  and  noise  analysis  is 
made  resulting  in  a  theoretical  maximum  performance 
criteria.   A  comparison  is  made  with  the  performance  of 
various  configurations  of  the  system  demodulator . 
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1.   Introduction., 

It  has  long  been  recognized  that  reliability  in  high 
frequency  radio  communications  is  limited  not  so  much  by 
the  power  which  can  be  radiated  by  an  antenna  as  it  is  by 
the  destructive  vectoral  addition  of  signals  arriving  at 
the  receiver.   The  signals  which  are  added  vectorially 
may  be  delayed  replicas  or  may  be  parts  of  two  different 
symbols.   In  both  cases  the  effect  is  the  result  of  the 
delays  in  multipath  propagation.   The  first  case  is  selec- 
tive fading,  and  the  second  case  is  intersymbol  interfer- 
ence. 

There  have  been  many  schemes  devised  in  attempts  to 
reduce  the  effects  of  multipath  propagation  phenomena. 
Several  forms  of  diversity  reception  including  time  civer- 
sity,  space  diversity  and  frequency  diversity  each  designed 
to  isolate  one  of  the  incoming  signals,  and  various 
combining  techniques  have  been  usedc  /l/  Other  schemes 
include  the  use  of  single  sideband  /z/   and  synchronous 
detection.  /3/ 

Probably  the  most  advanced  technique  in  combatting 
multipath  effects  is  the  "Rake"  system ,  /k/     This  system 
uses  a  cross  correlation  detection  scheme,  coupled  with  a 
method  for  isolating  the  signal  from  each  path,  and  then 
combines  the  signals  from  each  path  in  nearly  optimum 
proportions. 

All  these  schemes  have  their  own  advantages  and  dis- 
advantages.  Usually  the  best  results  are  obtained  by 
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those  schemes  requiring  large  antennas  or  complex  equip- 
ment,  or  large  land  areas.   Other  schemes  of  less  complex- 
ity can  expect  less  increase  in  reliability. 

The  modulation  scheme  analyzed  in  this  paper  is  one 
designed  to  operate  within  the  bandwidth  of  one  single- 
sideband  voice  channel  of  3  K0C0   Also  a  minimum  size  and 
weight  limitation  is  imposed  which  precludes  large  antennas 
or  space  diversity  techniques .   Instead  a  type  of  frequency 
diversity  is  employed  which  allows  the  transmitter  to  oper- 
ate at  a  nearly  constant  power  output .   A  coded  pulse  is 
used  for  both  the  "mark"  and  "space",  and  a  delay  line 
summer  is  used  to  add  coherently  the  output  pulses 0 

The  modulation  scheme  will  be  described  and  related 
to  the  principles  of  communication  theory  <>   A  signal  and 
noise  analysis  will  be  made  on  the  system  to  predict  error 
rates  for  various  signal  to  noise  power  ratios. 

The  methods  and  results  of  laboratory  testing  of  the 
demodulator  under  various  signal  to  noise  ratios  will  be 
presented „   Comparisons  will  be  made  between  the  predicted 
error  rates  and  the  actual  error  rates . 


2.   System  Operation. 

This  section  includes  a  description  of  the  basic 
modulation  scheme  and  the  method  used  for  demodulation 0 
The  principles  of  communication  theory  which  this  system 
implements  are  emphasized  where  employed . 

2.1  Basic  modulation  scheme. 

Quantized  Frequency  Modulation  (QFM)  is  the  name 
given  to  the  modulation  scheme  developed  by  the  Communi- 
cations Division  of  Hughes  Aircraft  Company.   This  scheme 
was  devised  in  an  effort  to  reduce  the  effects  of  multipath 
fading  on  the  reliability  of  teletype  or  mark-space  communis 
cations  systems. 

A  combination  of  several  techniques  is  applied ,  and 
these  techniques  will  be  described  in  later  sections. 

The  tapped  delay  line  approximation  to  the  matched 
filter  is  used  in  this  system  as  follows. 

A  combination  of  short  pulses  and  frequency  stepping 
is  used  to  reduce  the  effects  of  intersymbol  interference 
and  selective  fading  respectively.   The  short  pulses  are 
spaced  at  each  frequency  to  allow  time  for  the  path 
delays  to  die  out.   But  upon  detection  these  pulses  are 
added  coherently  by  means  of  the  tapped  delay  line.   Fig.  1 
is  an  idealized  graphical  representation  of  the  signal 
transmitted. 
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This  is  the  modulation  on  a  carriers   Notice  that  there 
is  uninterrupted  modulation  of  the  carrier,  and  hence  the 
output  power  is  constants   In  order  to  contain  the  modulation 
within  a  3KC  bandwidth  the  pulses  are  shaped  to  approximately 
a  cosine  squared  shape.   An  example  of  the  actual  shaped 
pulses  is  shown  in  Fig,  2  for  Mark— Space— Space—  Mark  as 
before . 

To  demodulate  this  Quantized  Frequency  Modulation 
a  demodulator  of  the  type  described  below  is  used,   Fig,  3 
is  a  block  diagram  which  will  be  used  in  the  following 
description. 

The  QFM  signal  is  obtained  from  a  standard  communica- 
tions receiver  and  is  put  through  the  four  filters,  but 
only  one  of  four  filters  will  pass  a  signal  at  any  instant 
of  time.  The  other  three  filters  will  pass  only  noise. 
The  amplifier  will  amplify  the  signal  to  a  maximum  level 
if  the  signal  is  large  enough  to  reach  that  level,  other- 
wise it  is  amplified  by  the  full  gain  of  the  amplifier. 
Each  channel  is  full  wave  detected  either  positively  or 
negatively,  and  the  channels  are  added  as  shown  in  Fig,  3» 
The  signals  are  then  summed  in  the  tapped  delay  line  and 
differential  amplifier. 

The  output  of  the  differential  amplifier  is  used  to 
establish  synchronization,  which  in  turn  gates  the  inte- 
grator to  operate  only  during  the  time  that  the  peak  signal 
occurs. 
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The  output  of  the  gated  integrator  will  be  positive 
or  negative  depending  on  whether  a  "mark"  or  "space"  code 
was  transmitted o   The  output  of  the  gated  integrator  is 
sampled  just  before  the  integrator  capacitor  is  discharged 0 
The  sampling  timing  is  also  controlled  by  the  synchroni- 
zation circuit o   The  polarity  of  the  integrator  output  will 
then  determine  whether  the  data  is  "mark"  or  "space"  and 

will  trigger  the  appropriate  data  flip-flop,, 

'  1 

Since  the  decision  is  based  solely  on  the  polarity  of 

integrator  output  the  decision  is  made  about  a  nearly 
zero  threshold o 

The  operation  of  the  system  has  been  described 0   The 
principles  of  communication  theory  employed  will  be  described 
where  applicable  in  the  following  sections » 
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2o2  Short  pulse  and  mult !•=  frequency  transmission 

The  effect  of  multipath  delay  for  pulses  which  are 
long  compared  to  the  delay  time  is  shown  in  FigQ  4o 

The  effect  of  multipath  delay  for  pulses  which  are 
short  compared  to  the  delay  time  is  shown  in  Fig0  5* 

It  can  be  seen  that  the  short  pulse  is  more  desirable, 
but  of  course  this  pulse  requires  more  bandwidth,  and  con- 
tains less  total  energyc 

Turning  now  to  the  selective  fading  effects  of  multi- 
path  propagation  and  directing  our-  attention  to  Fig,  6, 
it  may  be  noted  that  the  multipath  delay  is  different  for 
different  frequencies B   So  one  possible  approach  to  com- 
batting selective  fading  could  be  by  using  more  than  one 
frequency  in  the  transmission. 

In  this  QFM  scheme,  the  techniques  and  those  to  be 
described  in  Sections  2„3  and  204  are  combined.   In 
addition,  in  order  to  increase  the  total  energy  in  a  mark 
baud  and  space  baud,  ten  pulses  are  transmitted  for  each 
baud0   At  the  receiver  the  envelopes  of  these  pulses  are 
added  coherently  as  described  in  the  previous  section 0 
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2«,3  Coded  pulses  and  matched  filter  detection. 

It  has  been  shown  /5/  that  the  optimum  receiver  to 
detect  a  signal  in  white  Gaussian  noise  is  one  which 
calculates  the  likelihood  ratio  and  then  decides  if  a 
signal  was  present  or  not  depending  on  whether  the  like- 
lihood ratio  exceeded  a  given  threshold  value  or  not 
respectively o   The  receiver  which  does  this  is  a  cross 
correlator,  and  essentially  calculates  the  coefficient  of 
linear  correlation  between  the  received  signal  and  a  copy 
of  the  transmitted  signal 0   If  this  exceeds  a  specified 
value  then  the  decision  is  made  that  a  signal  was  present „ 
Since  this  is  a  continuous  process  Per) is  continuously  being 
calculated  as 


/_<  y^> 


TA 


C(r)    =-  T_79C   "F/  s,%)  j*#-v  *+       (o 


where  Si  (^   is  the  received  signal  and 

S^(t)      is  the  copy  of  the  transmitted  signal 0 
For  application  to  the  communications  problem  the 
decision  to  be  made  is,   "Was  the  signal  sent  a  mark  or 
a  space?" c   It  has  been  shown,  /l/  that  the  optimum 
decision  process  is  one  which  calculates  for  both 
the  mark  and  space  by  cross  correlating  the  incoming 
signal  with  the  "mark"  signal  as  transmitted  and  the 
"space"  signal  as  transmitted s  and  then  decides  which 
is  largerc 
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and  ft*  «**  (r;  -  £'-£0    T  Xr,    s<&  f*»*rK   £-T;  <(*       (*> 

where  Xr/««t  to  is  the  signal  representing  "Space" 
and  Sr*asM.(£>      is  the  signal  representing  "Mark". 

It  has  also  been  shown  /l/  that  the  probability  of 
making  a  wrong  decision  is  least  when   mark     space 
and  the  decision  threshold  is  zero.  The  cross  correlation 
of  the  "mark"  signal  as  transmitted  with  a  "space"  signal 
as  transmitted  for  this  case  would  be  -1. 

One  means  of  implementing  a  cross  correlation  detec- 
tor is  by  means  of  a  filter.  The  output  of  a  filter  for 
a  given  input  signal  is  the  convolution  of  the  signal 
with  the  impulse  response  of  the  filter. 


or 
to 


if  the  impulse  response  of  the  filter  h(t)  is  the  time 
reverse  of  the  input  signal,  X(tt),  then  the  output  is  the 
convolution  of  X(t)  with  X(-t).  Substituting  in  equation 
(5)  above 

CO 

#  ir)  *X  yi(P)  *  iy~ (~*j)  M         ^ 

P  C°  "f°°  *&  x  (t+v  ^  P 
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Since  \tf)   as  shown  in  equation  (1)  is  an  even  function 
for  periodic  signals 

few  =  f(-T;  (* ') 

(V,  --  f(-r,  .  &*.  ^y  v\t, ,  <*  ,t,  *        ''  '' 

In  any  physical  implementation,  integration  starts 

at  some  time  t0  and  is  carried  out  for  some  finite  time 

interval.   Hence 

C(T>    'J?  I?) 

for  any  physically  realizable  system. 

In  terms  of  probability,  the  question,  "Given  that 
a  signal  was  sent,  what  is  the  probility  that  it  was  a 
"mark"  and  what  is  the  probability  that  it  was  a  "space"  ?" 
must  be  answered.   Since  each  is  equally  likely,  the  a 
priori  probability  of  each  is  one-half. 

The  a  posteriori  probability  that  a  mark  (for  instance) 
was  sent  is,  from  /l/ 


tffp /S '*,*.*$)] fts*.** * £;]  +f>[U>/-Tsf^>]  F[Sjr«*(tJ J 


It  has  been  shown  /l/  that  for  mark  and  space  equally  likely 
the  decision  reduces  to: 
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"Decide   "mark"   if 

Jc     f(t)  S*a,K   ft)    d*  -Jo      f(t)  Syctcl  (t)  d*  L'Z> 

is  greater  than  zero  (positive).  Decide  "space"  if  (12) 
above  is  less  than  zero  (negative). 

If  a  truly  matched  filter  or  correlation  detector 
could  be  built  the  error  rate  has  been  shown  to  be 
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Figure  7 
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2.4  Tapped  delay  line  approximization  to  a  matched  filter, 

Turin  /6/  proposed  three  methods  of  implementing  a 
matched  filter.   One  of  these  was  the  tapped  delay  line0 
However,  he  pointed  out  that  there  are  physical  limitations 
to  this  approach.   The  filter  which  is  truly  matched  is 
matched  continuously  from  beginning  to  end  of  the  signal 
when  T=  0.   Fig.  8   shows  an  example  of  a  signal  and  the 
impulse  response  of  a  filter  matched  to  it. 


signal 


h(t)  for  filter  matched  to  signal 
Figure  8. 

A  signal  which  might  be  put  into  a  tapped  delay  line  and 
its  corresponding  matched  filter  is  shown  in  Fig.  9» 
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Figure  9 

This  particular  signal  pattern  chosen  is  one  of  the 
well  known  "Barker  codes"  /7,S/.   It  was  chosen  to  illus- 
trate a  desirable  property  which  is  possessed  by  some 
binary  codes,  namely  that  of  a  "good"  autocorrelation 
function.   Notice  that  the  interval  between  peaks  is  flat. 

However  if  instead  of  a  truly  matched  filter  a 
tapped  delay  line  is  used,  the  output  of  the  filter  is 
still  the  correlation  of  the  input  signal  with  the  impulse 
response  of  the  filter,  but  the  taps  are  not  as  "broad" 
as  the  pulse  length.   Instead  the  taps  could,  be  represented 
as  delta  functions  weighted  by  the  proper  sign  as  shown  in 
Fig.  10. 
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Figure  10 
Notice  that  the  output^  at  7  =  0  has  the  shape  of 
an  individual  signal  of  width  T  instead  of  the  triangle. 
The  delay  line  it  is  seen,  linearly  adds  the  voltages  at 
each  tap,  and  is  not  truly  "matched"  to  the  signal.   The 
next  section  will  deal  with  the  analysis  of  the  system, 
and  will  establish  a  mathematical  model  to  predict  theoret- 
ical  performance  maxima. 
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3.   Analysis  of  Demodulator. 

3.1  Mathematical  model  of  demodulator. 

In  order  to  analyze  the  performance  of  the  system 
and  to  be  able  to  predict  error  rates  for  various  signal 
to  noise  power  ratio  inputs  a  mathematical  model  was 
devised.  Models  have  been  established  and  performance 
analysis  done  for  several  of  the  well  known  binary 
communications  systems.  /9,10,ll/  However  most  of  these 
analyses  assume  a  signal  to  noise  ratio  of  3db  or  more0 
This  paper  deals  with  lower  signal  to  noise  power  ratios 
and  consequently  the  approximation  used  in  them  are  not 
valid  for  this  study.   The  signal  to  noise  power  ratios 
in  and  out  of  an  envelope  detector  have  been  described 
completely  by  Rice.  /12/  Techniques  in  the  application 
of  his  results  may  be  found  in  other  publications  as  well, 

A3,  H/ 

The  exact  expression  for  the  density  function  of  the 
output  of  an  envelope  detector  includes  a  modified  Bessel 
function.   At  this  point  the  high  speed  digital  computer 
was  used  to  obtain  the  signal  voltage  and  noise  power  for 
different  signal  to  noise  ratios.   This  system  is  essentially 
linear  from  input  to  the  full  wave  rectifier  of  each  channel, 
and  subsequent  to  that.   Hence  if  the  non-linearity  of  the 
full  wave  rectifier  could  be  analyzed  and  modeled  then  the 
entire  system  could  be  modeled.   The  amplifier  prior  to 
detection  will  be  assumed  linear,  and  section  two  demon- 
strated that  the  tapped  delay  output  was  the  linear  alge- 


braic  sum  of  the  signal  in  the  filter  at  every  instant 
of  time.   Hence  the  signal  to  noise  power  ratio  at  the 
input  to  the  full  wave  detector  will  be  taken  as  a  refer- 
ence.  Signal  to  noise  power  ratios  at  the  input  to  the 
demodulator  can  then  be  expressed  in  terms  of  the  refer- 
ence.  In  making  the  transition  from  full  wave  rectifier 
input  to  demodulator  input  it  will  be  assumed  that  the 
noise  figure  for  the  AGC  amplifier  is  zero, 

Another  assumption  in  the  model  is  that  the  signal 
pulse  is  a  square  pulse.   This  is  done  in  order  to  simplify 
the  analysis  of  the  model.   However  in  order  to  use  the 
model  to  relate  predicted  error  rates  to  actual  error 
rates  measured  a  conversion  factor  is  used.   Appendix  IV 
shows  in  detail  how  this  conversion  factor  is  calculated. 
Once  the  relationship  is  made  the  relative  positions  of 
the  predicted  curve  to  the  actual  measured  error  rate  curve 
is  established.   The  error  rate  versus  signal  to  noise 
ratio  per  pulse  at  the  input  to  the  detector  can  be  used 
as  a  reference,  and  error  rates  for  signal  to  noise  ratio 
at  any  point  before  the  detector  can  then  be  obtained.   This 
relationship  may  be  seen  by  comparing  the  graphs  of  Fig„ 
19  and  Fig.  20. 

The  model  is  considered  to  have  four  separate  channels „ 
Each  channel  has  an  ideal  filter  5$0  cycles  wide  centered 
at  6S0,  1260,  1940,  and  2620  cps.  The  signals  of  each 
channel  are  envelope  detected  and  added  together  as  shown 
in  Fig.  3.   The  output  of  the  tapped  delay  line  is  inte- 
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grated  for  one  half  the  pulse  period  during  the  time 
interval  that  the  summed  pulse  is  highest.   A  decision  is 
then  made  at  the  end  of  integration  time  on  the  basis  of  the 
polarity  of  the  integrator  output.   "Mark"  if  the  integra- 
tor output  is  positive  and  "Space"  if  the  integrator  output 
is  negative. 

The  probability  of  error  is  the  probability  that 
the  integrator  output  will  be  negative  (or  positive)  when 
a  "Mark"  (or  "Space")  signal  was  transmitted. 


Figure  11 
From  Fig.  11  it  may  be  seen  that  the  probability  of  error 
is : 


In  order  to  take  advantage  of  the  tabulations  of  the 
area  under  the  normal  curve  /15/  this  is  put  into  the  form 

r6  --  f  [>  -E^f  v] 

where  Erf  V  =     _J_         "  v 

rr    J 
V 


v; 


7.TT 


fVs*r-[  £:  it 
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3.2  Analysis  of  model. 

It  has  been  shown  by  Rice  /ll/  that  if  an  envelope 
detector  has  as  its  input  bandlimited  white  Gaussian  noise, 
the  probability  density  function  for  the  output  is  the 
Rayleigh  distribution.   Specifically  for  an  input  or 
bandlimited  white  Gaussian  noise  of  power  spectral  density 
NQ  watts  per  cycle  per  second,  the  output  of  an  envelope 
detector  has  as  its  probability  density  function  the  Ray- 
leigh distribution, 

f  (y)=   1    ye-y2/2  N2 

where  N2  ■  N0B. 

The  power  spectral  density  and  probability  density 
function  is  as  shown  in  Fig.  12.   As  a  signal  is  added  to 
the  noise  input,  the  density  function  shifts  to  the  right 
as  shown  in  Fig.  13  and  Fig.  14. 

The  probability  density  function  for  the  output  is 
now  a  function  of  the  signal  power  and  the  noise  power. 
Specifically  from  /14/ 

/•(y)=  JL  ts^c- <£.*,-   J.  If'^jr-) 

where  IQ  (x)  is  the  Modified  Bessel  function  of  the  first 
kind  and  zero  order  for  the  argument  x,  and  s2  is  the 
signal  to  noise  power  ratio  for  a  continuous  carrier. 
Most  treatments  of  envelope  detectors  will  make  a  closed 
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form  approximation  to  I0  (x)  for  large  or  small  signal  to 
noise  ratios,  because  the  density  function  becomes 

essentially  normal  for  large  signal  to  noise  ratios,  and 

2 

IQ  (x)  can  be  closely  approximated  by  e£  for  x<<1.   How- 

ever  for  this  study  these  approximations  were  crudest  in 
the  region  of  main  interest.   The  modified  Bessel  function 
I0  (x)  may  be  expressed  as  the  sum  of  an  infinite  series, 
namely 


<*>=  f 


»--°  £nwf 


Appendix  I  shows  the  infinite  series  solution  to  ob- 
taining the  mean,  and  variance  of  the  resulting  probability 
density  function. 

Appendix  II  is  the  computer  programs  used  to  calcu- 
late the  mean  and  variance  of  the  distribution  resulting 
from  different  signal  to  noise  power  ratios  at  the  input 
to  the  envelope  detector.   The  program  also  includes  some 
other  calculations  used  to  obtain  the  probability  of  error 
which  would  have  been  quite  laborious  if  done  by  hand. 
Also  included  in  Appendix  II  is  a  sample  of  the  program 
outputs. 

Fig.  16  illustrates  an  example  of  the  probability 
density  functions  for  all  four  channels  at  one  instant 
of  time,  and  how  they  are  combined. 

Although  the  distribution  of  the  sum  of  the  four 


22 


channels  is  not  truly  normal,  it  is  nearly  so.   After 
summing  ten  of  these  nearly  normal  distributions ,  by  the 
central  limit  theorem,  the  resulting  distribution  will  be 
very  close  indeed  to  being  a  normal  distribution. 

This  signal  is  then  integrated  for  an  interval  of  2 
milliseconds.   Appendix  III  describes  the  treatment  of 
integrating  a  signal  plus  noise.   For  a  continuous  carrier 
pulse  this  amounts  to  reducing  the  signal  to  noise  power 
ratio.   The  result  is  manifest  on  the  distribution  function 
as  a  change  in  the  ratio  of  the  D.  C.  signal  level  to  the 
standard  deviation  of  the  distribution,  and  the  probability 
of  error  is  reduced  accordingly.   Probability  of  error  is 
now 

PE  *  §     l-Erf  Q 


where  Q  = 


y  /  0(^7*- f  j  Tsj 


J^y   is  the  mean  of  distribution  of  the  output  voltage  of 
the  channel  containing  the  signal  during  one  pulse  dura- 
tion. 0\     is  the  variance  of  the  distribution  of  the 
output  of  the  channel  containing  the  signal  during  one 
pulse  duration,  or  fluctuation  noise  power „ 

cM  ?   is  the  mean  of  the  Rayleigh  distribution  resulting 

-  2 
from  envelope  detection  of  white  Gaussian  noise,  o   p  ^-s  t^1 

variance  of  the  Rayleigh  distribution  or  fluctuation  noise 

power  of  the  output  of  envelope  detection  of  white  Gaussian 

noise. 
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Probability  Density  Function  for  the  Output 
of  an  Envelope  Detector  for  Increasing 
Signal  to  Noise  Ratio.  From /1 1/. 

Figure   15 
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The  standard  deviation  of  the  Gaussian  noise  distri=> 
bution  was  taken  as  one,  and  thus  the  noise  power  (f^^   was 
one.   Then  the  signal  to  noise  power  ratio  is  just  the 
signal  power.   A  Control  Data  Corporation  1604  computer 
was  used  to  calculate  Q  and  the  probability  of  error 
obtained  from  tables  /l5/«   The  probability  of  error  for 
various  signal  to  noise  power  ratios  was  obtained  and 
plotted.   Curve  V  of  Fig.  19  is  the  curve  obtained. 

Laboratory  measurements  were  made  to  compare  actual 
error  rates  with  those  predicted,,   A  description  of  the 
apparatus  and  techniques  used  and  the  results  obtained  is 
contained  in  the  next  two  sections. 


27 


4.   Laboratory  Tests  of  Demodulator0 

A  laboratory  test  of  the  demodulator  was  desired  to 
compare  actual  performance  with  the  theoretical  predictions. 
In  order  to  do  this  it  was  necessary  to  devise  some  means 
of  detecting  and  counting  errors  and  counting  decisions 
made.   Appendix  V  is  a  description  of  the  error  detector 
which  was  built  up  from  Engineered  Electronics  Company 
digital  building  blocks.  Fig.  17  is  a  block  diagram  of  the 
laboratory  equipment  layout. 

The  automatic  gain  control  circuits  of  the  amplifier 
were  set  for  minimum  gain,  and  the  linear  range  of  the 
amplifier  was  determined.  Fig.  IS  is  an  example  of  the 
linear  range  of  the  amplifier. 

A  few  tests  showed  that  the  error  rate  was  quite 
good  for  good  signal  to  noise  ratios  but  seemed  to  have 
a  signal  to  noise  ratio  threshold  below  which  the  demodu- 
lator lost  synchronization.   Actual  errors  as  a  result 
of  wrong  decisions  were  still  quite  low0 

Investigation  of  the  synchronization  circuits  re- 
vealed that  for  the  signal  to  noise  ratio  at  which  the 
demodulator  lost  synchronization,  there  was  still  available 
a  reasonably  "clean"  signal  pulse  coming  from  the  differ- 
ential amplifier.   However  the  voltage  level  of  the  signal 
was  not  of  sufficient  amplitude  to  be  sensed  by  the  part 
of  the  synchronization  circuit  which  decided  if  a  signal 
(either  "mark"  or  "space")  was  present.  A  two  stage  ampli« 
fier  was  then  inserted  between  the  differential  amplifier 
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and  the  "signal  sensing"  portion  of  the  synchronization 
circuit.   Curve  II  of  Fig.  19  is  the  curve  of  error  rates 
measured  for  this  configuration . 

Visicorder  recordings  of  the  outputs  of  the  gated 
integrator  and  data  output  showed  that  additional  improve- 
ment could  be  obtained.   The  time  constant  of  the  operational 
amplifier  integrator  was  reduced 0   The  error  rates  obtained 
for  this  configuration  is  shown  as  curve  III  of  Fig0  19 ■ 

Vjsicorder  recordings  show  that  the  sensing  threshold 
was  considerably  above  zero.   The  input  resistance  to  the 
sensing  transistors  was  reduced  to  allow  the  required 
switching  current  to  be  reached  at  a  lower  threshold 
voltage.  This  was  the  final  configuration  tested,  and  the 
error  rates  for  it  are  shown  as  curve  IV  of  Fig.  19» 

In  Fig.  19  the  measured  signal  to  noise  ratio  was 
corrected  to  fit  the  model  as  described  in  Appendix  IVU 
Fig.  19  is  the  same  data  referring  the  model  signal  to 
noise  ratios  to  the  measured  signal  to  noise  ratio  per 
bit. 

As  discussed  in  section  3,  the  data  in  Fig„  19  and  20 
are  for  the  signal  to  noise  ratio  at  the  input  to  the  full 
wave  rectifier  as  a  reference,  in  order  to  refer  the  error 
rate  to  the  signal  to  noise  ratio  ner  bit  at  the  input  to 
the  filters  the  following  correction  is  required. 
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Each  filter  bandwidth  is  5^0  cps0   The  total  noise 
bandwidth  at  the  input  is  3  kc  -  400  cps  =  2500  cpse   There 
is  then  -^SS  or  4.5  times  as  much  noise  per  bit  at  the 
input  as  in  each  channel.   Signal  to  noise  ratio  must  then 
be  reduced  by  6.53db.   Fig.  21  shows  this« 

Any  measurements  made  in  the  future  on  this  system 
will  probably  be  done  by  a  true  reading  RMS  voltmeter  as 
was  done  in  our  experiments.   Figs.  19,  20,  and  21  are  on 
a  signal  to  noise  power  per  bit  basis .   Since  the  signal 
was  a  mark — space — mark — space,  etc,  the  average  signal 
power  was  just  one-half  of  the  power  per  bit.   Average 
signal  to  noise  power  ratio  is  then  3db  less  for  each 

rate  than  for  Fig.  21.   This  is  shown  in  Fig,  22, 
and  is  the  error  rate  for  the  signal  to  noise  ratio  which 
would  be  measured  at  the  input  to  the  demodulator,, 

Fig.  23  is  a  group  of  pictures  showing  the  QFM  signal, 
the  differential  amplifier,  output  for  two  signal  to  noise 
ratios  and  configurations. 
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5v/cm  .VFK  signal  (no  noise) 
differential  amplifier 
output 


•  5v/cn. 


10v/cm  liF'K  signal   +  noise 

S/li    (per  bit)  =  243db  (Ref  COS' 

1v/cm  output  of  differential 

amplifier 
PE  (II)  =  .16 


5v/cm  gated  integrator  output 
1v/ca  output  of  differential 

amplifier 
PE  =  .16(11) 


,5v/cm  differential  amplifier 

output 
S/K  (per  bit)  =  ,3db  (Ref  COS2) 
5v/cm  £;ated  integrator  output 
PE  (IV)  =  ,o.T9 


pV 


5v/cm  gated   integrator  output 
5v/cm  data   outout 
PE    (IV)   =  .OS! 


M  >  £ 
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5o   Conclusions. 

A  close  analysis  of  a  system  to  obtain  a  theoretical 
optimum  performance  level  provides  a  useful  standard  to 
which  actual  performance  may  be  compared 0   Such  an  anal- 
ysis may  show  the  possible  improvement  to  be  gained  by 
alterations  in  various  sections  of  a  system 0   Such  was 
the  case  for  the  QFM  system  studied 0   Curve  V  is  the 
theoretical  optimum  performance  criteria  for  this  system^ 
and  I  to  IV  are  performance  levels  for  various  configuations 
of  the  system.   It  appears  that  additional  improvements  in 
system  performance  are  still  possible .   However  these 
curves  cannot  show  the  limitations  of  the  synchronization 
circuit.   Perfect  synchronization  and  zero  threshold  are 
assumed. 

In  any  physically  realizable  system  synchronization 
is  not  perfect  and  the  threshold  is  some  small  value  above 
and  below  zero. 
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APPENDIX  I 

Infinite  Series  Solution  to  Mean  and  Variance  of  the 
Probability  Density  Function  for  the  Output  of  an 
Envelope  Detector  for  Various  Signal  to  Noise 

Poxtfer  Ratios 

From  /ll/  the  probability  density  function  for  the 
output  of  an  envelope  detector  is 

where 

S2  is  the  signal  to  noise  power  ratio  at  the  input 
to  the  envelope  detector , 
N2  is  the  noise  power  at  the  input  and 
I0(Y)  is  the  modified  Bessel  function  of  the  first  kind 
and  zeroth  order  of  the  argument  Y  . 
I  (Y)  can  be  expressed  as  the  sum  of  an  infinite  series 

The  mean  of  this  distribution  is 


<« 


-co 


since  the  distribution  has  only  positive  values  for  X  and 


c?  *  (A 

Term  by  term  integration  yields 


41 


*y 


For  digital  Computer  Operation  this  series  may  be 
nested  as 

The  computer  program  which  does  this  is  in  Appendix  II, 

The  Second  Moment  for  the  distribution  function  has 
been  evaluated  in  closed  form  as 

The  variance  or  fluctuation  noise  power  is  thens 

(fup   and  J(ma.y   be  evaluated  for  different  signal  to  noise 
power  ratios. 
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Appendix  II 
The  digital  computer  program  written  to  compute  the 
mean  and  variance  of  the  distribution  for  the  output  of 
envelope  detector  is  shown  below.  Sample  outputs  are  also 
included  in  this  appendix. 


PROGRAM  SIGNOIS 
;      SIGNAL  VOLTAGE  AND  NOISE  POWER  FROM  AN  ENVELOPE  DETECTOR 
DIMENSION  TRM( 100), EMI (90), EM2( 90), XOISPR( 90), SI GPR< 90) 
i»££5I£  90),EM1SQO(90),TDCSIG(90),TNOSPR(90),ERF<90) 
2,DB(90),SIGINT(90),ERFINT(90) 
.  P0WER=0 

DO  20  1=1,90 
'  IF(  1-10)4,4,5 

4  P0WER=P0WER*.01 
GO  TO  10 

5  IF(I-40)6,7,7 

6  P0WER=P0WER+.1 
GO  TO  10 

7  IF(I-60)8,9,9 

8  P0WER=P0WER+.5 
GO  TO  10 

9  POWER=POWER+1.0 

10  SIGPRU  )=P0WER 
XEGPR=-SIGPR(I) 
T=1.0 
SUM=1.0 

TRM(1  )=1.0 
DO  15  N=2,100 
T=T+1.0 
„   0TRM(N)=((TRM(N-1)  )  •  (  S  IGPR(  I )  »  (  2.0»T-1  .0)  )  )/(2.0«(  (  T-l  •  <»•)(  T-1  .0)  ) 
15   I SUM=SUM+TRM(N)  ' 

EMI (I )=SUM»EXPF(X EGPR)« 1,2533 141 
EM1SQD(I)=EM1(I)tEM1 (I) 
EM2U  )=2.0»(1.0«-SIGPRU)  ) 
XOISPRt I)=EM2U)-EM1SQDU) 
DCSIG(I)=EM1( I )-l. 25331 41 
TDCSIG( I)=10.0»0:SIG( I) 
TN0SPRU)=10.0»UOISPRU)*1.28761102) 
QPR=TNOSPR( I)  s 

20  ERFU  )=TDCSIGU)/SQRTF(QPRK 
-  PRINT  21 

21  OFORMAT  (97H   SIG  TO  NOISE  IN     FIRST  MOMENT  OUT     SECOND  MOMEN 
l?OUT    DC  SIGNAL  VOLTAGE    NOISE  POWER  OUT) 

PRINT  3,  (SIGPRU)  ,  EMI  (I )  ,  EM2  (  I  ) ,  DCSIG  (I )  ,  XOISPR  (I )  , 
1  1  =  1,90) 
3     FORMAT  (5F20.8) 
PRINT  29 
PRINT  98,  (SIGPRd  ),TDCSIGU),TNOSPRU),ERF(  I)  ,  I»l,  90) 

98  FORMAT  (4F20.8) 
PRINT  29 

29    FORMAT(IHI) 

DO  30   1=1,90 

DB(  I)=10.0*  LOGIOF(SIGPRU)) 

SIGINTU)=1.2»T0:SIGU) 

30  ERFINTU)*  SIGINT  (I )  /SQRTF(  TNOSPRU  ) ) 
PRINT  31,  (DBU),SIGINT<I),ERFINTm. 1-1,90) 

31  FORMAT  (3F20.8) 
END 

END 

PRINT  99, (TRM(I), 1-1,100) 

99  FORMAT  (F20.8) 
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APPENDIX  III 
FILTER  CHARACTERISTICS  OF  AN  INTEGRATOR 
The  impulse  response  of  an  ideal  integrator  which 
integrates  for  a  time  interval  T  is  shown  as  a  solid  line 
in  Fig,  A0  /17/ 


r 


Impulse  response  of  an  ideal  integrator. 
Figure  III-A 

To  obtain  the  transfer  characteristics  in  the  frequency 
demain  translate  the  coordinates  an  amount  T/2,  shown  as 
dotted  in  Fig.  A.   Transform  this  to  the  frequency  domain,, 


H  (v; 
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A  delay  in  the  time  domain  of  T/2  seconds  is  manifest 
as  a  phase  shift  in  the  frequency  domain  of  <f *  „   Hence 


The  power  output  of  such  an  integrator  would  be 
Po     ~  /        \Hii^)  \     C-{lj)    4-f 

where  G(W)  is  the  power  spectral  density  of  the  input 
signal o 

The  power  output  of  the  integrator  becomes 

ft./  \^^Xl0       u 


~0 

CO 


^g-^^  *"  Tf  {  T        A  ,  , 

rrfr       (r{^J)    *{ 


for  one  sided  frequency  spectrum0 

For  this  problem  T  ■  2  x  10  -*   seconds  „   The  frequency 
of  the  first  zero  is  then 


k& 


For  any  one  of  those  three  channels  which  for  any- 
instant  contain  noise  only  the  power  spectral  density  is 
as  shown  in  Fig.  13  from  /ll/ 

From  the  Rayleigh  distribution,  the  variance  or  fluctuation 
noise  power  is  .429/N^. 

Then      .^^"-  £  Izc  **-»<>) 


There  will  always  be  three  channels  of  noise  power 
,429tfij2  an(i  the  noise  power  from  the  signal  channel  going 
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into  the  integrator  which  will  add  directly  since  the 
noise  from  different  parts  of  the  spectrum  are  not  coherent, 
The  noise  power  in  the  signal  channel  will  be  increasing 
with  the  increasing  signal  power „   Thus  the  noise  power 
may  be  considered  as  four  channels  of  noise  pov/er  plus 
the  increased  noise  power  in  the  sign£l  channel  due  to 
the  signal  times  noise  products.   This  increase  may  be 
closely  approximated  by  curve  fitting  techniques  applied 
to  Fig0  III-B.   By  linear  regression  the  increase  in 
noise  power  is 

where  X  is  the  signal  to  noise  power  ratio0 

This  increase  in  noise  power  is  essentially  flat  over 
half  the  original  bandwidth  of  the  noise  input  to  the 
detector e   Thus  the  Noise  Pov/er  Spectral  Density  into  the 
integrator  will  be 

and  the  pov/er  output  of  the  integrator  is 
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The  power  spectral  density  input  to  the  integrator 
will  be  as  shown  in  Fig„  III-C. 
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Figure  III-C 

In  order  to  obtain  a  tractable  model  the  noise  power 
spectral  density  will  be  assumed  to  be  flat  from  zero  to 
B/2  cycles  per  second  with  a  total  noise  power  equal  to  the 
original  noise  power.   The  noise  power  out  of  the  detector 
then  is 

t  }*(»){  ii 


p.  - 


-co 
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Transforming  variables 

""A 


-*?> 


-  .  7  7  V  /V 


This  modifies  the  error  rate  to  be 
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APPENDIX  IV 

The  analysis  of  the  QFM  demodulator  assumed  a 
square  pulse  shape  for  each  bit  transmitted  as  shown  in 
Fig.  lo   The  pulse  transmitted  is  shown  in  Fig0  2o   In 
order  to  bridge  the  gap  between  the  model  and  the  actual 
modulator  the  following  conversion  proceedure  was  applied 0 

The  power  in  a  constant  unmodulated  carrier  of 
amplitude  A  is  A2/2. 


The  cosine  squared  pulse  which  is  the  output  of  the 
rectifier  is  obtained  from  an  input  to  the  rectifier  of 
a  waveform  like  a  carrier  100  per  cent  modulated  by  a 
sine  waveG   The  power  in  this  signal  is  3/2  ~     <>   In 


terms  of  the  peak  voltage  of  the  modulated  carrier 

2     * 


^.--J-.i  fa 


5"  5  (r)   -  s 


The  "mark"  or  "space"  decision  is  based  on  the  polarity 
of  the  integrator  output  after  2  millisiconds  of  integration 
The  output  of  the  integrator  for  the  modulated  input  is 


/  _ 


% 


The  amplitude  of  a  constant  carrier  pulse  for  an  equivalent 
voltage  output  of  the  intigrator  would  be  a&2   PK. 
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The  power  in  this  carrier  is 

The  power  in  a  constant  carrier  signal  which  will  produce 
the  same  output  at  the  integrator  as  a  modulated  carrier 
with  a  power  per  bit  of  Pg  is 

p         p        r.  6?*^  /V" 

Since  all  signal  to  noise  ratios  are  referred  to  a  noise 
power  of  one,  the  Signal  Power  is  the  signal  to  noise 
ratio,  and 

I  0     L  OQ.     Ppq    *    /0   Loa    \.?Z  +  10  i-tp  f& 


Thus  2«53db  should  be  added  to  the  signal  power  per  bit 
as  measured. 

To  refer  the  model  to  the  measured  signal  power  2»53db 
must  be  subtracted  from  the  theoretical  noise  power., 
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APPENDIX  V 
ERROR  DETECTION  INSTRUMENTATION 
Some  means  of  detecting  errors  was  necessary  in 
order  to  obtain  an  error  count 0   Since  the  output  data  was 
digital,  a  digital  system  was  used  for  error  detectionD 
Engineered  Electronics  building  blocks  were  arranged  as 
shown  in  the  diagram  below. 


RECEIVED 
DATA 


DEMODULATOR 
TIMING 


_F.F#B   | 


Counter 


Flip  flop  A  assumed  the  state  corresponding  to  the  "Mark" 
or  "Space"  data  received.  Flip  flop  B  assumed  the  state  of 
the  transmitted  data.   The  outputs  were  gated  according  to 
the  following  equation, 

C .■  (A  .  B)   -   (A  .  B) 

Error  =  (C  .  D) 
D  was  derived  from  the  interface  timing  circuits 0 
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